Power Counting and Perturbative One Pion Exchange in Heavy Meson Molecules 
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We discuss the possible power counting schemes that can be applied in the effective field theory 
description of heavy meson molecules, such as the X{3872) or the recently discovered Zi,(10610) and 
Zi, (10650) states. We argue that the effect of coupled channels is suppressed by at least two orders 
in the effective field theory expansion, meaning that they can be safely ignored at lowest order. 
The role of the one pion exchange potential between the heavy mesons, and in particular the tensor 
force, is also analyzed. By using techniques developed in atomic physics for handling power-law 
singular potentials, which have been also successfully employed in nuclear physics, we determine 
the range of center-of-mass momenta for which the tensor piece of the one pion exchange potential 
is perturbative. In this momentum range, the one pion exchange potential can be considered a 
subleading order correction, leaving at lowest order a very simple effective field theory consisting 
only on contact-range interactions. 

PACS numbers: 13.75.Lb,14.40.Lb,14.40.Nd,14.40Pq,14.40Rt 
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I. INTRODUCTION 

Heavy meson molecules are a long established theoreti- 
cal prediction of hadronic physics [ll-ly] • The discovery of 
the X(3872) by the Belle collaboration 0, and the subse- 
quent confirmation by CDF DO and BABAR 
has provided so far the strongest candidate for a bound 
state of heavy mesons. Owing to the closeness of the 
X(3872) to the D°*D° threshold, an interpretation in 
terms of a bound state of these two mesons is both ap- 
pealing and natural. In turn, the recent finding of the 
^6(10610) and ^6(10650) resonances by Belle [B, [12, 
which are located just a few MeV above the B*B and 
B*B* thresholds respectively, also calls for a molecular 
description. Moreover, from heavy quark spin symme- 
try O, [U ■^fi should expect other low- lying states, the 
partners of the Z^'s, in the bottom sector [l6{. 

The shallow nature of the aforementioned candidates 
for molecular states, apart from being instrumental in 
their identification, also indicates the existence of a sep- 
aration of scales between long and short range dynamics. 
The heavy mesons are far apart from each other and con- 
sequently are not able to resolve the details of the short 
range interaction that may be ultimately responsible for 
binding them. We therefore expect that heavy meson 
molecules will be amenable to an effective field th eory 
(EFT) description of their properties and decays [l7l - [l9{ . 

In the EFT formalism, the long range interaction be- 
tween two heavy mesons (composed of a heavy and a 
light quark) is constrained by the low-energy symmetries 
of the system, in particular chiral symmetry, which de- 
termines the pion-exchange dynamics. In turn, the short 
range interaction is mimicked by local, contact-range op- 
erators or counterterms. A second ingredient of the EFT 



description is the existence of a power counting, that is, 
an organizational principle from which operators can be 
ordered from more to less relevant. From power counting 
we expect to be able to express operators and observables 
as a power series in term of a small expansion parameter 
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where vq is the order at which the EFT expansion starts, 
I'max the order at which we perform the calculation, and 
xq, the expansion parameter, can be written as the ratio 
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with Q (Ao) the generic energy scale associated with the 
long (short) range physics. A priori power counting en- 
tails a remarkable advantage: the relative error of a cal- 
culation of order v is known to be Xq^^ . Of course, the 
realization of this promise of accuracy depends on the 
implementation details and on the unambiguous identifi- 
cation of the expansion parameter. This manuscript will 
try to deal with this problem. 

At this point a natural question arises: how important 
are the pion exchanges in the description of the heavy 
meson molecules? Of course, having a bound state re- 
quires the non-perturbative treatment of a certain sub- 
set of the meson-antimeson interaction, but not neces- 
sarily the one-pion exchange (OPE) potential. A nice 
illustration is provided by the X(3872), in which the dis- 
tance between the heavy meson and antimeson is so large 
{\/{r^ ~ 10 fm) that even pions may not be clearly dis- 
tinguished. This means that we can use a contact-range, 
pionless EFT in which pions are already considered to be 
short-ranged ^20,, 21} ■ In such a case we are left with a 
much simpler theory. Even though the hard scale is the 
pion mass, Aq = m^, which is an extremely light scale in 
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hadronic physics, the soft scale Q can be even lighter if 
the molecular state is close enough to threshold. In par- 
ticular, Q can be identified with the wave number of the 
bound state, that is, Q — y/2fij^B, where fj,ji is the re- 
duced mass of the two particle system and B the binding 
energy. For the D*^D^ mesons conforming the X(3872), 
Q is of the order of mere tens of MeV. Thus the EFT 
expansion is expected to converge fairly well, providing 
a motivation and explaining the success of contact-range 
descriptions of the X(3872) ^Hiil- 

However, the applicability of pionless EFT is subjected 
to limitations. In the X(3872) we can list two: on the 
one hand, the pion exchanged between the D*° and 
meson is almost on the mass shell, and consequently its 
effects spread over a much larger range than expected. 
That is, the irrelevance of pion exchanges is not so evi- 
dent as it appears at first sight. Luckily, as shown in X- 
EFT [l9j, pious are actually perturbative in the X(3872). 
On the other, if we consider the charged D*'^D~ com- 
ponent of the X{3872), which is essential for explain- 
ing certain branching ratios [1^, , the associated soft 
scale is Q ~ 125 MeV, of the order of TOtt. Thus the 
explicit inclusion of the charged channel may lie outside 
the range of applicability of a pionless EFT. If we now 
consider the bottom sector, the Zb(10610) and Zf,(10650) 
resonances are a only few MeV away from the B*B and 
B*B* threshold Q. However, the large reduced mass of 
these systems imply that Q ^ corresponds to a bind- 
ing energy of merely 4 MeV at which the wave function 
will start to probe the pions, requiring a non-perturbative 
theory jH]. 

In all the previous examples the inclusion of pion ex- 
changes may be required for a proper EFT description 
of the molecular states. This can be rather cumbersome, 
owing to the rich angular momentum coupled channel 
structure triggered by tensor forces, especially in the 
D*D* and B* B* cases. However, there exists a binding 
energy range in which the OPE potential is subleading 
with respect to the contact range interactions and hence 
perturbative. In this energy window the pionfull EFT 
will still be a contact-range theory at the lowest (or lead- 
ing) order (LO), as pion exchanges will not enter until 
next-to-leading order (NLO). In the two-nucleon system 
the paradigmatic example of this kind of EFT is the Ka- 
plan, Savage and Wise (KSW) counting [H, Is^] , which 
served as inspiration for X-EFT [l9|. 

The question we want to answer is: where does the 
limit between perturbative and non-perturbative pions 
stand? If we take a second look to the deuteron in the two 
nucleon system, which shares many similarities with the 



heavy meson molecules, we see that the convergence of 
the EFT with perturbative pions is numerically marginal. 
As shown by Fleming, Stewart and Mehen (EMS) [3l| 
by a thorough next-to-next-to-leading order (N^LO) cal- 
culation, the convergence of the KSW counting in the 
deuteron case is limited at most to Ao ^ 100 MeV, a 
rather low figure (indeed smaller than m^). Taking 
into account that in the deuteron the wave number is 
7 = 45 MeV, the previous breakdown scale translates 
into a rather slow convergence rate. In fact, the theory 
would not converge at all had the deuteron be bound by 
about 4 MeV or more. The purpose of this work is there- 
fore to find the corresponding Smax below which pions 
are perturbative, as this limit will gives us essential in- 
formation about the convergence of the EFT for heavy 
meson molecules. 

Of course, performing a N^LO calculation is beyond 
the scope of the present paper. Apart from that, the 
scarce experimental input available on heavy meson- 
antimeson systems makes the previous task impractical: 
we do not have the information required to fix the coun- 
terterms. Therefore, we need to resort to a more indirect 
path. In the two-nucleon system the solution was pro- 
vided by Birse [s^l. The idea is to consider the OPE in 
the chiral limit, in which it reduces to a pure power-law 
potential of the type The interesting thing here 

is that these kind of potentials have been studied and 
analyzed in detail in the field of atomic physics. The 
long-range solutions of OPE in the chiral limit can be 
expressed in terms of a a particular type of series, as 
shown by Cavagnero [s^ and Gao [s^] • Birse was able to 
show that one can extract the breakdown scale of a the- 
ory containing perturbative pions with the help of these 
techniques, and the results help to explain remarkably 
weU the lengthier N^LO calculations of EMS t31i]- In this 
work we will extend the observations of Birse to the pe- 
culiarities of the tensor force in heavy meson molecules. 

The article is structure as follows: in Sect.lUwe present 
a brief overview of certain key EFT ideas that will help 
to put the results of this work into proper context. In 
Sect, mil we will write the OPE and contact range po- 
tentials between a heavy meson and antimeson, and in 
Sect. IIVI we will discuss the different power counting 
schemes that we can apply. In Sect. |V] we will extract 
the breakdown scale of perturbative pion theories, and 
finally, in Sect. |Vl]we will discuss the results and their 
implications on the EFT treatment of heavy molecular 
states. We have also included an Appendix containing 
the technical details involved in the derivation of the EFT 
potential at lowest order. 



It is interesting to notice that even though the original experi- 
mental analysis by Belle [Til . IT^ locates the two states above the 
B* B and B* B* thresholds, this depends on the Breit-Wi gne r 
parametrization employed for the Z^'s as stressed in Ref. |27ll . 
This work also indicates that the two Z(,'s may be located be- 
low threshold and hence be bound states (instead of two-particle 
resonances) after all. 



II. GENERAL CONSIDERATIONS 

The purpose of this section is to provide a quick re- 
view of the EFT formalism for non-relativistic two-body 
systems. The discussion is heavily based on the EFT for- 
mulation of the two-nucleon system (see Refs. jS 5143 91] for 
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reviews), which can be trivially translated and applied 
to heavy meson molecules with minor modifications, as 
demonstrated in Ref. [28]. Thus, we begin with naive 
dimensional power counting, as originally proposed by 
Weinberg [io, IHl, and then explain the modifications 
that have been required to successfully formulate a non- 
relativistic EFT to two- body swtems forming shallow 
bound states [13, [H HO, SI 13 ■ As we are dealing 
with heavy meson molecules, we will include notation 
related to heavy quark symmetry. In what follows toq is 
the mass of the heavy quark conforming a heavy meson, 
and H — P, P* is used to denote a generic heavy meson 
with orbital angular momentum I = between the heavy 
quark and the light quark. The heavy meson P (P*) has 
total spin s — {s — 1) and hence it is a pseudoscalar 
(vector) meson. If we are specifically dealing with the 
charm or bottom sector we will particularize the P, P* 
heavy meson notation by D, D* and B, B*. We also use 
D(*) and B(*) as a generic for the pseudoscalar/vector 
cases. We will only consider the case in which the light 
quark is the u or d quark. The extension to the strange 
sector is straightforward. 



A. Power Counting 

The formulation of EFT depends on the existence of a 
separation of scales: we can distinguish between Q, the 
low energy scale that characterizes the physics we are 
interested in, and Ag, the high energy scale at which the 
effective description we are using stops to be applicable. 
In the EFT framework, the two-body potential can be 
expanded as a power series on the ratio of these scales, 
leading to 



^^EFT - J2 ^'"^ + ^ 
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where vq > —1 is the order at which the expansion begins 
and fmax the order at which we perform the calculation. 
In chiral (nuclear [13, El and heavy hadron Q ) EFT 
the generic scale Q usually includes the momenta of the 
two interacting particles and the mass of the pion, while 
Aq refers to the mass of the rho meson or the momentum 
scale at which the internal structure of the two particles 
starts to be resolved. 

As can be seen, we expect the theoretical error in the 
determination of the EFT potential to decrease, at least 
as long as Q < Aq. At this point we should take into 
account that the light scales include on the one hand 
the pion mass ttItt, which does not change 0, and on the 
other the momenta Q ^ p,p' of the two heavy mesons, 
which can vary. This means that the expansion of the 
potential (and the scattering observables) is only valid for 



sufficiently small p.,p'. To avoid the related divergences 
with taking p, p' > Aq in loops we usually include a cut- 
off A that serves as an intermediate scale between Q and 
Aq (i.e. we take Q < A< Aq). 



B. The Scaling of Operators 

The power counting assignment of a certain operator 
tells us about its scaling properties. A contribution to 
the (momentum space) potential is assigned the order v 
if it scales as 



y('')(Ag) = A'^y('')(Q) 



(4) 



under the rescaling of all the light scales by a factor A. 
This means that if we reduce the momentum or the pion 
mass by a factor i? = j, with i? < 1, thus increasing the 
separation of scales by a factor of A = > 1, the size of 
the order v contributions will decrease as R"^ . That is, 
such a contribution becomes smaller the larger the scale 
separation, justifying their power counting assignment. 

Scaling is very interesting in the sense that it deter- 
mines the behaviour of the potential in coordinate space. 
If we consider the order v contribution potential in mo- 
mentum space, we simply have[f| 



FM(A<r,Ag) = A''FM(q,Q) 



(5) 



where we have now explicitly considered the dependence 
on the momentum exchanges between the two particles, 
q = p — p' . After Fourier-transforming into coordinate 
space, the previous scaling translates to 



A3+"^y('')(r,g), 



(6) 



which admits two kind of general solutions, contact range 
and finite range. The contact range solution is trivial to 
construct from the Dirac (5-function and its derivatives, 
yielding 



V^''\r,Q) = C,d''5{f). 



(7) 



where d denotes a general derivative of the Dirac 6. Of 
course, parity constraints imply that this kind of contri- 
bution only appears at v = 2n. On the other hand, the 
finite range solution must comply to the form 



F('')(r,Q) 



(8) 



with f{x) an arbitrary (non exclusively power- law) func- 
tion, which decays exponentially at large distances as it 
stems from meson exchanges. 



Unless we are considering chiral extrapolations. 



^ Notice that we have restricted ourselves to the local potential 
case. Non-localities only appear at high orders in the chiral ex- 
pansion (t' = 4 in the two-nucleon case) and in addition they are 
suppressed by the mass of the heavy mesons. 
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In the previous equations Ci, and are constants with 



dimensions of 1/ [energy] . From power counting we 
expect the related energy scale to be Aq, that is, 

1 . „ 1 



A, 



u+2 



and F„ 



^^0 



(9) 



However, as we will see in the following paragraphs, 
sometimes there is an unexpected contamination of Ci, 
and Fi, by a light scale. This will trigger a change in 
the power counting of the potential, promoting a certain 
contribution from a higher to a lower order. An illustra- 
tion is given by the scaling arguments discussed in the 
introduction of Ref. [4a . 



C. Iteration 

If we are considering scattering states, we should plug 
the potential into the Lippmann-Schwinger equation, 



T^V + VGoT , 



(10) 



where V is the EFT potential, T is the T-matrix and 
Go = 1/{E — Hq) the resolvent operator. Analogously, 
if we are considering heavy meson molecules we should 
iterate the potential in the bound state equation 



(11) 



to obtain the binding energy and the wave function |5'_b). 

Within the EFT framework the previous equations 
are expected to be re-expanded according to the power 
counting of the potential. In this way, we are guaranteed 
to be able to estimate the error properly. For this, we 
need to take into account that the scaling of the resol- 
vent operator is given by 



Gn 



Mi?' 



(12) 



where fin is the reduced mass of the two body system, 
which for a non-relativistic theory can be considered a 
high energy scale, ~ Aq. From the point of view 
of power counting, the scaling of Gq means that only 
contributions to the EFT potential such that 0{GoV) ^ 
1 should be iterated, i.e. the order piece of the 

potential. 

The existence of shallow bound states between heavy 
mesons imply that there is a contribution to the poten- 
tial of order Q^^. However, when we compute the EFT 
potential from the chiral lagrangian we only find contri- 
butions of order and higher. The OPE potential is 
order Q", as can be trivially checked from 



(13) 



as has been discussed extensively in nuclear EFT [20|, |29|, 
Isol . . A way to see this is to consider the coupling 

Go of the lowest order contact operator: if there is a 
bound state, Go cannot be perturbative, that is, small. 
We can account for the fact that Go is big by multiplying 
the naive expectation for its size by a big number, say 
Ao/Q, yielding 



Go(Q) 



QAo 



(14) 



from which it is obvious that Co{Q) is of order Q^^, 
promoting the lowest order contact range potential from 
order Q*^ to Q~^. A different way is to solve the bound 
state equation with a contact operator, in which case we 
obtain that Go should scale as 



Go(Q) 



1 



7A0 



(15) 



with 7 = ^/2fj,ii B the wave number of the heavy meson 
molecule, and B the binding energy. For a shallow bound 
state we have j Q and the wave number is the light 
scale contaminating Go. 

Moreover, there is a second mechanism that can lead 
to the promotion of the potential, in particular OPE, 
to order Q~^. It is based on a well-known argument 
of Ericson and Karl j.6] about the relative strength of 
OPE in heavy meson molecules. The idea is that the 
intrinsic strength of the OPE potential, i.e. Fq as defined 
in Eq. (|9]), is roughly independent on the heavy quark 
massQ As a consequence, the ratio 



V^(o)Gol/(") m Q 



Ao Ao 



(16) 



is only small for the naive identification /i/j ^ Aq, for 
which the ratio above scales as expected if the OPE po- 
tential is a Q° operator. However, while the chiral hard 
scale Aq is fixed, the reduced mass of the two heavy me- 
son systems scales as mg in the heavy quark limit. Even- 
tually we can have jJa S> Ao, enhancing the pion loops 
to the point of making the previous ratio to be of or- 
der Q". In such a case the OPE potential will become 
non-perturbative, as expected from Ref. @. 



D. Coupled Channels 

If we consider heavy meson molecules, we expect from 
heavy quark symmetry that the energy gap between the 
thresholds of the three different combinations of pseu- 
doscalar and vector mesons, that is 



HH = PP,PP*,P*P* 



(17) 



and it is also clear that contact range potentials must be 
at least of order Q'^. This means that naive dimensional 
analysis is not enough to determine accurately the order 
of an operator if there is a bound state. What is missing 
is the unexpected contribution from a light energy scale. 



^ Actually, it is proportional to the axial pion coupling g^, which 
contains a contribution of order plus corrections of order 

rrig^, rriQ^ and so on. 
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will shrink as Ag ~ ^/mq [4^, |43|- This indicates, at 
first sight, the necessity of a coupled channel approach. 
In such a case, the Lippmann-Schwinger equation can be 
written as 



Tab = Vab + Vac Gq Tcb , 



(18) 



where A,B,C = PP,PP*,P*P*. If we are considering 
the scattering in the vicinity of a given threshold a, then 
we set A = B = a in the Lippmann-Schwinger equation. 
In addition, we set the center-of-mass energy to zero E = 
at the a threshold. This means that we must include 
a proper energy shift for the resolvent operator if Gq is 
acting on a channel C a, that is 

G^\E) \q,C) = {E-£-- Ag) I q, C) , (19) 

where we have considered the inverse of the Gq operator 
for simplicity. In the expression above, /ic is the reduced 
mass of channel C and the energy shift is given by Ag, — 
Ma — Mc, with Ma and Mc the total mass of channels 
a and C respectively. 

However, the counting of the Go operator is not nec- 
essarily Q on a channel G ^ a. If we consider the ratio 
of the Go operator evaluated at the C and a channels 
respectively, we find that 



G^AE) 



2^^A E~q^ 



Ag) 



Q 
Ac 



^ (20) 



where Q ~ yJ2^BE ^ q (we are assuming similar reduced 
masses), and Ac — y^2/icAg is the momentum scale 
related to the coupled channel effects. We can distinguish 
here two possible situations: (i) Ac ^ Ao and (ii) Ac ~ 
Aq. If Ac <C Ao, we can consider the coupled channel 
scale to be light, Ac ~ Q, and the resolvent operator 
will be of order Q, as expected (a more formal account 
can be consulted in Refs. [H, H^). On the contrary, if 
Ac ^ Ao, then the resolvent operator for the channel 
G ^ a will be suppressed by two powers of Q/Ac, that 
is, by two orders in the EFT expansion. 

In the particular case of heavy mesons, Ac = 
•,y2/ZHHAQ, where /zhh is the reduced mass of the sys- 
tem and Aq is the energy split between different pairs 
of heavy meson systems. For the charm sector, Ac — 
520 MeV for the DD - DD* and DD* - D*D* pairs and 
Ac ^ 740 MeV for DD - D*D*. For the bottom sec- 
tor, the coupled-channel scales are very similar, yielding 
Ac ^ 490 MeV and Ac ^ 700 MeV respectively. The 
previous estimations indicate that Ac is similar to the 
hard scale of the theory, Aq ~ 0.5 — 1.0 GeV, which indi- 
cates that coupled channel effects are suppressed. 

A surprising aspect about the coupled channel effects 
is that the associated momentum scale is similar in the 
charm and bottom sectors. This can be easily under- 
stood if we consider the heavy quark limit, mg — >■ oo, in 
which we expect ^hh ~ and Aq ^ l/iriQ, meaning 
that A 



c 



That is, the suppression of the cou- 
pled channel effects is basically independent of the heavy 



quark mass. Even though in the heavy quark limit the 
heavy mesons are degenerate, the energy split between 
them remains beyond the scope of a non-relativistic chi- 
ral EFT. The convergence of the EFT description for 
non-relativistic particles depends on momentum scales 
instead of energy scales. However, from the scaling of the 
center-of-mass energy with respect to the heavy quark 
mass, i.e. E ~ p^/niQ, we see that the energy window 
in which the EFT is valid decreases as the heavy quark 
mass grows. 



III. THE LEADING ORDER POTENTIAL 

The lowest order potential in the effective field theory 
description of heavy meson molecules can be decomposed 
into a contact-range and a finite-range piece 



^ HH 



(0) 



(21) 



where the subscripts c and f are used to denote the con- 
tact and finite-range character of each contribution to the 
potential. In the equation above we have labelled the po- 
tential by the expected naive dimensional scaling, rather 
than the actual scaling that depends upon the particular 
power counting under consideration. Of course, the exis- 
tence of a shallow molecular scale will imply that at least 
Vc is of order Q~^. However, we will not consider the 
actual order of the different contributions of the potential 
(and their consequences) until the next section. 



A. The Contact Range Potential 

At LO we will naively expect a total of six countert- 
erms in a given isospin channel. This number correspond 
to the number of s-wave channels for the HH system: a 
0++ state for PP, two opposite C-parity states, 1+"'" and 
1+- for the_P*P/PP* system, and the three 0++, 1+" 
and 2++ p*p* states. However, the contact range is con- 
strained by the requirements of heavy quark spin sym- 
metry (HQSS) [H, [13, which in turn implies that the 
six counterterm figure is reduced to only two indepen- 
dent counterterms [3]. We will call these counterterms 
Goa and Gob. Ignoring the coupled channel effects, the 
potential reads (see the Appendix for details) 



Goa , 



(0) 



C,P'P/PP 



(0) 



. {q^ 1^ ) — Goa — Gob , 

.((f,l++) -Goa + Goh, 



C,P*P/PP 
Vl?l,p,{q,Q++)^Goa-2Gob: 

V^(5°^.p.(g,2++)=Goa + Go6, 



(22) 

(23) 
(24) 

(25) 
(26) 
(27) 
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depending on the J^*^ quantum number and the specific 
HH system under consideration. The expUcit represen- 
tation of these two counterterms in the cou pled channel 
basis of HH states can be consulted in Ref . [16| . 



B. The Finite Range Potential 

The OPE potential between a heavy meson and anti- 
meson is local, and hence it only depends on the ex- 
changed momentum between the heavy mesons: 



(28) 



with q = p — p' . The explicit consideration of the pseu- 
doscalar / vector heavy meson channels allows us to write 
the potential as 



1^F?P^PP(9)=0, 



(0) 



F,P'P^pp*(9) — 2/2^-^'^^ «2 ^ ,,2 



(0) 



F,P*P*->P*P 



q + ^l' 
^ ^ Si - qS2- q 



(29) 
(30) 

, (31) 



depending on the particular heavy meson channel un- 
der consideration. In the equations above g is the pion 
axial coupling, — 132 MeV the pion decay constant, 
mjr = 138 MeV the pion mass, ti(2) the isospin oper- 
ators on the heavy meson 1(2), £1(2) is the heavy vec- 
tor meson polarization, and 5*1(2) is the spin operator 
for intrinsic spin S — I. For the P*P — >• PP* chan- 
nel, we use /i instead of m^r: as a consequence of the 
different masses of the P and P* mesons, the pion is 
emitted with the zero component of the momentum dif- 
ferent to zero. In the static limit this is equivalent to 
changing the effective mass of the pion. We have that 
^2 — ^2 _ .^j^j^ ^Yie mass splitting between the 
P and P* mesons. In the heavy quark limit, mq — > 00, 
we obtain /i = ra^ -I- ©(I/ttiq), and it would be practical 
to expand the potential in powers of 1/toq. While this 
is the situation we find in the bottom sector (BB*), the 
charm one (the DD* potential) is more complicated as we 
have iJ? < 0. The OPE potential thus acquires a small 
imaginary piece that is related to the probability loss in- 
duced by the open decay channel DD*/DD* — >• DDtt. In 
this regard it is not clear up to what extend the static 
approximation to the potential holds. The rigorous the- 
oretical treatment of this situation requires the inclusion 
of pions as dynamical degrees of freedom and the explicit 
consideration of the DDtt three body channel |50|] , which 
indicates a significant impact of the three body dynamics 
on the X(3872) DDtt decay rate, but only a mild ef- 
fect on the DD* wave functions (that is, the residue of the 
X(3872) pole in the language of Ref. [SO])- This seems to 
indicate that, if we are only interested in the wave func- 
tions, we can simply ignore the imaginary piece of the 
OPE potential for fi^ < and continue using the static 
limit. 



Apart from the previous, there are additional pieces of 
the potential that mix the different heavy meson chan- 
nels. However, as we have seen, coupled channel effects 
can be safely disregarded (even in the heavy quark limit) 
as their size is similar to the short range effects beyond 
chiral symmetry. Therefore we have decided to ignore 
the particle coupled-channel terms. 

The previous potential has been computed for the 
meson-meson case. For obtaining the meson-antimeson 
potential in the isospin symmetric basis, we perform a G- 
parity transformation: this changes the sign of the P*P* 
potential, but leaves the p*p/pp* potential unchanged. 
However, the p*p/pp* potential is better written in a 
definite C-parity basis, for which we employ 

\P*Pir^))^l= [\P*P)~^\PP*)] , (32) 

where rj is the intrinsic C-parity of the meson-antimeson 
system. This translates into a factor of —77 for the 
P*P/PP* potential when expressed in this basis. The 
final form of the meson-antimeson potentials is thus 



(0) 



F,P*P(?7) 



5 - ei-qe2 -q , . 
^^ = ^^/^--- ,2 + ^2 ^ (33) 



(0) g ^ ^ 

p^p.p.io) ^ -^2 



Si ■ qS2- q 

q^ + ml 



(34) 



The general form of the potential can be schematically 
written as 
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Qi ■qb2-q 

q^ + yP' 



(35) 



where oi, &2 represent the particular spin operators and 
/i the particular value of the effective pion mass to be 
used in each case. In this form, the intrinsic C-parity 
is then to be taken rj — 1 for the P*P* channel. 

The coordinate space potential can be obtained by 
Fourier-transforming the momentum space one, in which 
case we obtain 



'^f!hh(^) = n • T2 ai • 62 



riT2-Tl 



di ■ b2Wc{r) 



+ {3 di ■ f b2 ■ f — ai ■ b2) Wrir) 



(36) 



where Wc (r) and Wt (r) are the central and tensor 
pieces of the potential, which read 



Weir) = 



247r/2 fir 



Wrir) 



2, ,3 



.9 M 



247r/2 jir 



1 



3 

jir 



(37) 
(38) 



The interesting feature of the tensor force is that it can 
mix channels with different angular momentum. 



Or, more properly, the spin-spin part of the potential. 
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C. The Partial Wave Decomposition 



with 



In this subsection we will show the explicit partial wave 
decomposition of the potential. Wc will use the spectro- 
scopic notation '^^'^^Lj for characterizing a certain par- 
tial wave with spin 5, orbital angular momentum L and 
total angular momentum J. If two or more partial waves 
are coupled we will indicate it with a dash, for example 
5P2-5F2 or ^Fs-^Ps-^Fs-^Hs. 

There are certain rules on how partial waves coupled 
to each other. On the one hand, the central operator, 
which we define as 



a 



p.p. 
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(r) 



€1-62 , 

Si ■ S2 , 



(39) 
(40) 



conserves parity, C-parity, spin, orbital and total angular 
momentum. Consequently it does not mix partial waves. 
On the other hand, the tensor operator, given by 



(f) = 3ei-re2* 



5f;P*(f) = 3 5l•f52•f■ 



el • £2 ! 
■ Si ■ S2 , 



(41) 
(42) 



conserves parity, C-parity and total angular momentum, 
but neither spin nor orbital angular momentum. How- 
ever, parity and C-parity conservation imply that the 
tensor force can only change the spin S and orbital an- 
gular momcntiim L by an even number of units. This 
in turn implies that we have four possibilities, which we 
define as follows 



'J. 
'{J 

'J. 



J 



'(J + i)j, 

- 2)j - "Jj - 



(43) 
(44) 

(45) 
(46) 



where {U)C stands for (un)coupled. In this notation, ^Ci 
is ^Si-^Di, 5(72 stands for ^P2-^F2 and ^-^Cg would be 
^F^ ^Pt, ^F^ ^H^. Of course, low J coupled channels 
can be uncoupled or have less members than expected. 
Examples are the ^Pq (^C'o) or the ^Sq-^Dq (^-^Cq) par- 
tial waves. 

The matrix elements of the central and tensor opera- 
tors between different ^'^+^Z/j partial wave are schemat- 
ically calculated as 



(ff {fJ,b\0i2{f)\fJ,a) 



S.L.J S A 



o 



.hJa"MtMa ^.J{SbSaLbL^) > 



(47) 



where o, b are the initial and final partial wave, charac- 
terized by the vector \{Sa(b)La(b))Ja(b)Ma(b))-, and {/z} is 
whatever internal angular momentum quantum numbers 
a and b have. The two Kronecker (5's are a consequence 
of total angular momentum conservation. In the case of 
the central operator, which also conserves S and L, we 
simply have 



(48) 



C 



pp' 

J{1L) 



= 1. 



^J{SL) — 2 ['^('^' " 



1)-41 



(49) 

(50) 



in the PP* and P*P* cases respectively. 

We write the matrix elements of the tensor operator di- 
rectly in one of the four coupled channel basis previously 
defined. For the uncoupled channel case we have 



sr\'uj) = 



-1, 
+1. 



In the ^Cj case we obtain 



s^f^ecj) 



-2, 

SZfCCj) = +2, 



(51) 
(52) 



(53) 
(54) 



1 



2J+1 

J-1 



-3^J{J + l) 
1 



2J+1 
J -I 



-3^J{J + 1) 



-3VJ(J+1) 
J + 2 



-3v/J(J + l) 
J + 2 



(55) 



(56) 



The other two set of coupled channels, ^Cj and ^~^Cj, 
only happen for the P*P* system. In the ^Cj case, we 
obtain the matrices 



szrec) = -i, 



(57) 



'j>2 



2J+1 

J + 5 3a/^(^ + 1) - 2 
3^yj{J + l)-2 J-4 



(58) 



where there is no J = partial wave with fits into the 
^C.j scheme (the ^Pq wave is unphysical). In the ^~^Cj 
case, we obtain 



-V2 
-V2 -2 



2 

75 V 5 

2 _ 7 Vl 

5 5 

6 VI _8 I 
5 5 5 / 



(59) 
(60) 
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C) 



3J(./-1) 
(2J+1)(2,/-1) 



2J(J+1) 



(2J+3)(2J-1) 



/ 3(J+l)(./+2) 
\ V (2J+3)(2J+1) 



3J(.7-1) 
(2J+1)(2J-1) 

2.7-4 
2,7-1 



6(J-l)(J+l)(2J+3) 
(2J+1)(2J-1)^ 





2,7(J+1) 



(2J+3)(2J-1) 
6(J+l)(.7-l)(2.7+3) 
(2.7+l)(2.7-l)^ 
(2J+5)(2J-3) 
(2J+3)(2J-1) 
6J(J+2)(2,/-l) 
(2J+3)^(2.7+l) 



3(J+l)(.7+2) \ 
(2J+3)(2J+1) ' 



6J(J+2)(2J-1) 
(2J+3)^(2J+1) 
2J+6 

2J+3 , 



For J = 0, 1 the number of angular momentum channels 
is smaller than expected and the corresponding tensor 
matrices are written in the ^Sq-^Dq basis for J — and 
iPi-^Pi-^Fi for J = 1. 



IV. THE POWER COUNTING MAP FOR 
HEAVY MESON MOLECULES 

The presence of a bound state between a heavy me- 
son and antimeson indicates the necessity of promoting a 
piece of the EFT potential from order to Q^^. If the 
molecular state is shallow, as happens in the X{3872) 
or the Zf,(10610) and Z;,(10650), probably it is enough 
to promote a contact range interaction only. However, 
if the state is sufficiently deep (exactly how deep will be 
the subject of discussion of the next section), we are re- 
quired to additionally promote the OPE potential as well. 
In this second case, the renormalizability of the theory 
will in turn generate important changes in the countert- 
erm structure already at LO, as has been repeatedly dis- 
cussed in the EFT description of nuclear forces [s^, IsTI - 
[ssj . Therefore we distinguish two cases: power counting 
with perturbative and non-perturbative pious. 



A. Counting with Perturbative Pions 

We begin by considering a theory in which there is at 
least a Q^^ contact operator but where pions are pertur- 
bative. In this case, pions represent a NLO correction. 
As there are two contact operators at LO, Cqo and Cob, 
we have three options in what pertains to operator pro- 
motions, defining three power counting schemes: 

(a) Coa is of order Q^^, 

(b) Cob is of order Q^^, 

(c) both Coa and Cob are of order Q~^. 

Power countings (a), (b) and (c) are the EFT restatement 
of the observation made by Voloshin J^] about the heavy 
quark spin structure of the Zi, resonances: depending on 
which piece of the heavy quark spin symmetric interac- 
tion is responsible of the appearance of the two low-lying 
Zh states, we should expect a total of four or six molec- 
ular s-wave states of B(*^B(*^. The later situation, six 



states, correspond to the promotion (a), while the for- 
mer, giving four states, with (b) and (c). It is worth 
commenting that the third possibility (c) is quite general 
and can also accommodate other situations. For example 
if Coa < 0, Cob > and Coo + Cob is negative and attrac- 
tive enough as to bind, then there will be six states, but 
at most only two of them will be shallow (the 1++ and 
2++). However, this specific situation does not neces- 
sarily correspond to what we encounter in the Zb(10610) 
and Zb(10650) resonances. 

Power countings (a) and (b/c) also differ with respect 
to the treatment of coupled channel effects. The contact 
range interaction responsible for mixing different particle 
channels is Cob. In contrast, Coa is diagonal in the HH 
space. This means that in the (b/c) counting the cou- 
pled channels effects are promoted by one order in the 
EFT expansion, and are therefore 0{Q^), instead of the 
naive estimate 0{Q^). In turn, the coupled channel ef- 
fects dependent on the OPE interaction will still be of 
order Q^. 

In what regards the subleading contact operators, C2m 
their scaling is affected by the corresponding scaling of 
Co. By this we mean that if Co is (Q°)5 then C2n 
is (<3^"), as can be deduced from renormalization 

group analysis arguments [i^. Of course, only the C2„ 
operators with a similar HQSS structure as their 
counterparts will be promoted, depending on the whether 
we are considering the (a), (b) or (c) countings. However, 
the point is that, independently of the kind of promotion, 
we will have new free parameters at NLO. This also 
implies that without new data to fix the C2's we cannot 
proceed to NLO. Finally, the expected relative accuracy 
of a LO calculation is of order Q/Aq. 

There is still one contribution that we have not ad- 
dressed: the effect of heavy quarkonia states with the ad- 
equate quantum numbers in the vicinity of any of the HH 
thresholds. This situation may happen in the X{3872) if 
the Xci{2S) lies nearby, as has been discussed in Ref. 56']. 
However, the coupling to heavy quarkonia is surely going 
to represent a subleading contribution only. The nat- 
ural expectation is that the counterterm mixing these 
two states is of order Q'^. Therefore the coupled chan- 



Strictly speaking in (b) we should only expect three low-lying 
states, the two 1^ states and the p*p* 0+"*" state. However, 
if we recover the coupled channel effects, a fourth PP 0+"*" state 
may be expected. In this regard, it is more natural to use the 
promotion (c) to obtain the four states. 
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nel effects will be of order Q*^ and hence subleading if the 
thresholds of the X(3872) and the Xcii'^S) are very close, 
or order otherwise. This type of contribution, being 
subleading, does not really affect the LO calculation and 
its accuracy. 



B. Counting with Non-Perturbative Pions 

As we will see in the next section, pion exchanges 
can eventually become non-perturbative above a cer- 
tain critical energy in heavy meson- antimeson systems. 
Interestingly, the non-perturbative treatment of pion 
exchange will be able to significantly alter the power 
coun ting scheme, as happens in the two-nucleon sys- 
tem [llllll-lii]. We notice that Ref. [H] already devel- 
ops an EFT description of heavy meson molecules with 
non-perturbative OPE. However, the previous work only 
considers the particular case of the p*p/pp* channels. 
In the next lines we will concentrate instead on the gen- 
eral HQSS structures relating the different HH particle 
channels. 

We must begin with distinguishing which piece of the 
OPE potential is to be iterated, the central or the tensor. 
If only the central pieces is non-perturbative, the changes 
in the power counting are minimal with respect to the 
theories in which we already iterate the contact range 
operators, as demonstrated by Barford and Birse (57| . 
The reason is that the 1/r singularity of central OPE 
is not enough to change the power-law behaviour of the 
wave function at short distances, which is intimately re- 
lated to the scaling of the counterterms (see, for example, 
the detailed discussion of Ref. 53]). In such a case, we 
simply refer to the results of the previous section. 

On the contrary, the non-perturbative treatment of 
tensor OPE entails significant changes in the scaling of 
the contact range operators, as has been extensively dis- 
cussed in the context of nuclear EFT [H [HHH, Hi [11] . 
The reason is that tensor OPE is a singular potential 
behaving as l/r^ at distances below the pion Compton 
wavelength. The most evident effect of a singular poten- 
tial is on the LO counterterms: if the singular potential 
is attractive, the renormalizability of the theory requires 
the inclusion of a contact interaction to stabilize the cut- 
off dependence [H, This can be appreciated if we 
consider the solution of the reduced Schrodinger equation 
for a — potential, which for short enough distances 
read 



Uk{r) 













0,3 




(as. 


1 sm 










. r 



(62) 



The problem with the wave function above is that it is 
always regular at the origin, no matter what the val- 
ues of the semiclassical phase is. If we try to regular- 
ize the potential by cutting it off at short distances, i.e. 
V{r; Tc) = V{r) 9{r — Vc), we do obtain a value of ip{rc)- 
However, there is no well-defined — limit for (^(r^): 
the semiclassical phase simply oscillates faster and faster 



on its way to the origin. The solution is to make the 
counterterm oscillate 16(1, 



A^flC'o(?-c) 
27rr? 



3/4 

— I — I cot 



2 

as 



— + ip 
r 



(63) 



so ip remains a constant. In contrast, if the potential is 
repulsive, we have a unique and regular solution at short 
distances 



— exp 

as J 



03 

r 



(64) 



and no counterterm is required [58|, l59| . 

The (angular momentum) coupled channel case can be 
a bit more complicated however. In this case the impor- 
tant factor are the eigenvalues of the potential matrix 
in coupled channel space: if there are n attractive (i.e. 
negative) eigenvalues, we will need n{n + l)/2 countert- 
erms [58, 59]. For the tensor matrices we have calculated 
in the previous section, the eigenvalues are 



App(0++) = 0, 



App.(r 

App.(l^ 



r{2, 



-2} 
-1} 



Ap.p.(0^ 
Ap.p.(l^ 
Ap.p.(2^ 



n/3}, 



t{1 + V3,1 
r{2,-l}, 
T {1,-2, 1 + 73,1 



(65) 

(66) 
(67) 

(68) 
(69) 

V3} , (70) 



with T — Ti ■ T2, that is, we have included the — 77T1 ■ T2 
factor multiplying the potential of Eq. (1551) for obtain- 
ing the overall sign right. From the eigenvalues above it 
is apparent that the Coa and C06 counterterms are not 
enough as to renormalize the scattering amplitude if ten- 
sor OPE is non-perturbative. 

Each distinct negative eigenvalue requires a different, 
independent counterterm. However, if two different chan- 
nels share a negative eigenvalue, a common counterterm 
will be able to renormalize the two channels (60j . In prin- 
ciple we can count a total of eight counterterms for renor- 
malizing the six s-wave states (one per channel, with the 
exception of the 2'*"+ that requires three counterterms). 
This figure is then reduced by counting the shared eigen- 
values, that is, (i) the 1+" P*P and P*P* molecules can 
be renormalized with the same counterterm, (ii) the 0"''^ 
and 2++ P*P* share another eigenvalue and finally (iii) 
another common contact operator between the 1++ and 
2"'""'" channels. In total we end up with five independent 
counterterms. 

All this makes the theory with non-perturbative pions 
rather cumbersome. In addition, with five contact in- 
teractions HQSS, although still there, is not so manifest 
as in the pionless theory 0. It is interesting to notice 



What is actually happening is that certain s- to d-wave operators 
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though that the l+~ P*P and P*P* states are still ex- 
pected to be degenerate: this situation corresponds to 
the quantum numbers of the 2'f,(lQ610) and Zf,(10650) 
resonances. Luckily, as we will see in the next section, 
non-perturbative tensor forces are only expected in the 
isoscalar B(*)b(*) states. 

However the non-perturbative OPE counting has cer- 
tain advantages too. The momentum dependent C2 op- 
erators are not expected until order Q^l"^ at least (or even 
Q^: this is still an unresolved issue in nuclear EFT [ssf). 
The lowest subleading order corrections are expected to 
happen at order Q^: they are subleading contributions 
to the Co operators that are needed for absorbing the ex- 
pected divergences in (particle) coupled channels, owing 
to the Xjr^ singularity of the tensor force. Thus the ex- 
pected relative error of a LO calculation is now (Q/Aq)^, 
representing an improvement over the (Q/Aq) error in 
the perturbative pion case. 



V. THE PERTURBATIVE TREATMENT OF 
ONE PION EXCHANGE 



is expected to scale as 



Ac' 



(71) 



indicating the breakdown scale of a theory with central 
perturbative pions, where it is understood that we only 
consider the central piece of y(o). The evaluation is triv- 
ial for s- waves and p = 0, for which the generic low en- 
ergy scale Q can only be identified with the pion mass 
TOtt, leading to [3l| 

1 24^/2 



Ar 



where cr = 5 S{S + 1) - 2 and t = 2 /(/ + 1) - 3 with 
S and / the total spin and isospin of the system. In 
general, this number is quite large for isovector channels 
(Ac 1 GeV) , so the central piece of the OPE potential 
can always be treated as a perturbation in these cases. In 
contrast, for the 5' = isoscalar case, corresponding to 
a I^{JP^) = 0+(0++) P*P* molecule, the central piece 
can be quite important, but not so much as the tensor 
piece, as we will see later. 



(72) 



A. Central One Pion Exchange 



B. Tensor One Pion Exchange 




ia) 



FIG. 1. Feynman diagrams corresponding to (a) the OPE 
potential between two heavy mesons and (b) the first itera- 
tion of the OPE potential. The ratio of diagram (b) over (a) 
can be used to determine the energy range in which OPE is 
perturbative. 

The perturbative character of the central piece of the 
OPE potential can be determined by a direct comparison 
of the diagrams of Fig. ([TJ. This corresponds to com- 
puting the ratio of the matrix elements of the operators 
y(o)Gol/(o) and V^°\ which in terms of power counting 



are being promoted from to Q~^. However this may require 
the additional promotion of all (s- to s-wave) counterterms with 
two derivatives. Moreover, non-perturbative pions can change 
the bound state spectrum considerably from the expectations 
based on the HQSS of the potential, see Ref. [gij for an explicit 
example. A further issue, which we do not discuss in this work, is 
the possibility to go beyond the additivity of the long and short 
range forces, see footnote 6 of Ref. [6^ . 



The determination of the breakdown scale for a power 
counting in which tensor OPE is order and hence 
perturbative is not trivial. The reason is that the ar- 
gument of comparing the relative size of diagrams does 
not work when considering the tensor piece of OPE. In 
first place, tensor OPE does not directly act on s-waves, 
but only indirectly owing to transitions to intermediate 
d-wave states. This means that we should consider sec- 
ond and third order perturbation theory to obtain the 
breakdown scale 



(0), 



Ai 



(73) 



How- 
in the 



where by V^^^ we refer to the tensor piece only, 
ever, there is a serious problem at this point: 
equation above the numerator and denominator are di- 
vergent. The reason lies in the divergent behaviour 
of the tensor force at distances below the Compton wave- 
length of the pion, which induces a linear (quadratic) di- 
vergence in the denominator (numerator), see for exam- 
ple Ref. [iSl for the details of this kind of calculations. 
From the power counting point of view, V^^'^GoV'^^^ re- 
quires the inclusion of the order Q correction to the Cq 
counterterm, while V'^'^^GoV^-^^GoV'-'^^ involves in addi- 
tion two new order counterterms, a purely s-wave one 
and another one connecting the s- and d-wave channels 
for the result to be properly rcnormalized 0. The diver- 
gence spoils the evaluation of the previous ratio. If we 



* Of course, this is assuming that the leading counterterm appears 
at order Q°. If we are considering instead a two heavy meson 
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include counterterms, the ratio is again finite, but is also 
contaminated by whatever we assume to be the hard scale 
of the counterterms. In principle it is still possible to find 
the value of a posteriori by means of a full calcula- 
tion as in the nuclear case 31], but that requires plenty of 
experimental information about the low energy scatter- 
ing of heavy mesons to fix the counterterms. Therefore 
we need to resort to other kind of arguments. 

The failure of standard perturbation theory for in- 
verse power-law potentials (1/r") is well-known in atomic 
physics, where techniques have been developed to handle 
this type of potentials. Of particular interest in this re- 
gard is the work of Cavagnero [s^l , who pointed out that 
the divergences of the perturbative series are analogous 
to the role of secwZar perturbations in classical mechanics, 
that is, small perturbations that however end up diverg- 
ing at large enough time scales. The solution is to formu- 
late a secular series in which we redefine (or renormalize) 
some quantity, in this case the angular momentum, in 
order to obtain finite results at arbitrary order. By this 
we specifically mean that the zeroth order perturbative 
approximation to the Z-wave radial wave function is taken 
to be 

(74) 

instead of the customary solution r~^/^ J;_^i (/cr), where 
v(ka) is the renormalized angular momentum, k the 
center-of-mass momentum, a is the coupling constant and 
Jn{x) the Bessel function of order n. The interesting fea- 
ture is that at low energies the secular expansion can be 
reinterpreted as a particular resummation of the pertur- 
bative expansion, i.e. the renormalized angular momen- 
tum v{ka) can be expanded as a power series in a. This 
observation, if translated to the EFT potential between 
heavy mesons, indicates the existence of an energy range 
for which the inverse power-law potential can be treated 
perturbatively at long distances Q 

For the particular inverse power-law potential of inter- 
est in this work, the 1 /r'^ potential (corresponding to the 
form of the tensor force for m^r < 1), the secular series 
for the wave function have been studied and analyzed in 
detail by Gao [34] for the uncoupled channel case. These 
techniques were extended by Birse [32| to the coupled 
channel case, with a particular emphasis on the applica- 
tions to nucleon-nucleon scattering. In the present work, 
we particularize the results of Birse [s^l to the specific 
type of coupled channels appearing in two heavy meson 
systems. 



system with a shallow bound state, then the leading counterterm 
is order and the C2/C4 of order Q^/Q^. 

In this statement we are implicitly assuming the presence of ex- 
ponential suppression (e"""") for standard perturbation theory 
to make sense at long distances and of suitable counterterms to 
absorb the divergences we commented previously. 



In agreement with EFT expectations, the secular ex- 
pansion of the wave function converges fast at long dis- 
tances / low energies. Meanwhile, at high enough ener- 
gies a very interesting thing happens: the secular expan- 
sion becomes a non-converging series in powers of the 
coupling constant above a certain critical value of the 
center-of-mass momentum 0. That is, the secular ex- 
pansion can no longer be considered as a perturbative 
series in disguise. From the EFT viewpoint this critical 
momentum is to be identified with the hard scale of the 
theory: if we map the secular series of the wave function 
at low energies onto the EFT expansion 

*HH^E^*^^^(|:)^ (75) 

it is apparent that the mapping will stop being valid at 
Q — At, where Q is to be understood as the momentum 
(as the pion mass is fixed). This means that At is the 
tensor breakdown scale we were looking for. Above this 
critical value of the momentum the EFT description of 
the two heavy meson system will require the full iteration 
of the tensor force. 

In the following we will consider the OPE potential in 
the /i — > limit, for which it simplifies to a pure power- 
law potential of the type \/r'^: 

2MHHV(°)(r) = ^S,, (76) 

where Sj is the tensor force matrix, 03 is a length scale 
related to the strength of the tensor force, and the sub- 
script jjH denotes any combination of heavy meson and 
antimeson for which the OPE potential is non-zero. As 
can be appreciated, only the tensor piece of the OPE po- 
tential survives in this limit. Even though /i = can 
be physically interpreted as taking the chiral and heavy 
quark limit simultaneously (m^r — > and mg — )■ cx)), here 
we simply consider the /i — > limit as a mathematical 
convenience that allows us to directly apply the results 
of Refs. [H, 113] • Deviations from — will in general 
rise the value of the critical momentum, as they will re- 
duce the strength of the potential at distances fir > 1 
and therefore increase the perturbative character of the 
interaction at long distances. Analogously, the finite size 
of the heavy mesons is expected to weaken the poten- 
tial at short distances. In this regard, the fj, — (plus 
Aq = 00) determination are merely a lower bound: the 
real critical momenta will be higher. We will estimate 
the finite pion mass corrections to the critical momenta 
at the end of this section. 

Curiously there is a case for which the condition fi ~ 
is almost fulfilled: the A'(3872), in which m^r and Aq 
are very similar. In such a case the critical momenta 
calculated in this work will really represent the actual 
boundary to the perturbative treatment of tensor OPE. 



Specifically i^ix) encounters a non-analyticity in a; = fe a. 
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1. The Uncoupled Channel Case 

We start with the uncoupled channel case, for which 
the tensor operator Sj = Sj is simply a number (instead 
of a matrix in the coupled channel basis). In this case, 
the reduced Schrodinger equation reads 



— u 



k.J 



5. 



«3 , Ijilj + 1) 



(77) 



where Ukj is the reduced wave function, Ij is the angu- 
lar momentum of the two body system, which for the 
uncoupled channel case coincides with the total angular 
momentum of the system {Ij = j), and k is the center- 
of-mass momentum. We have |S'j| — 1 in all the HH 
uncoupled channels, with the exception of the ^Pq par- 
tial wave for which \Sj\ = 2. 

The general solution of this equation (for k ^ oEB) can 
be written as the linear combination 



Uk,] (r) = a (r; k) + (3 rjj (r; k) , 



(78) 



where ^k,j [f) and rjkj (r) are expressible as series in terms 
of Bessel functions 



4-00 

E 



+ 00 



Ki VrJ,n+^{kr) , 



(79) 



{r-k)^ (-1)™ W^J-™-.(fcr). (80) 



m— — oo 



The shift z/ is a function of the dimensionless parameter 
K = ka^. The bm coefficients can be expressed in terms of 
a recurrence relation that admits an explicit solution [s^ j . 

The dependence of the renormalized angular momen- 
tum v in terms of n is determined by the condition that 
i> is the zero of a characteristic function 

Fi{iy,K)=0, (81) 

where the function Fi(iy, k) is in turn defined by 

Fiiiy, k) = (i.2 _ ^1) - 'f. _ Roi-,,)] , (82) 

where = (^j + 5)- In the equation above, i?o(j^) can 
be expressed again in terms of a recursive relation 



(i/ + l)[(l.+ l)2-I.2]_^2 «2^^^^(^) 



(83) 



Birse [32| notices that between 20 and 30 iterations are 
enough as to obtain v^n) with six significant digits. As 



For simplicity we are expressing ui^ j in terms of wave functions 
which do not have a well-defined fc — > limit. This limitation can 
be easily overcome by adding energy-dependent normalization 
coefficients, see Ref. |34|| for the details of such a normalization. 



already commented, the analytical properties of the func- 
tion i'(k) are the essential ingredient for determining the 
range of applicability of the perturbative reexpansion of 
the wave function defined in Eqs. ^E^, 1^ and ([50]) . 

The behaviour of the solutions ly = can be de- 

scribed as follows: for k ~ (plus the condition A; 7^ 0) we 
have the Schrodinger equation for free waves with angu- 
lar momentum I. In such a case, the solution of Eq. (|77|) 
are the Spherical Bessel functions, which is equivalent to 
having 



j^(k = 0) = i^o = + 



1 



(84) 



If K is smaller than a critical value Kc, the shift can be 
written as 



i/(k < Kc) — 1^0 ^ S v . 



(85) 



with < S ly < ^. For sufficiently small k the 5 v cor- 
rection is of order \k\ (k^) for Ij — {Ij 7^ 0), where the 
exact expression can be consulted in Ref. [32] . For larger 
values of k, 5v is still a power series in k and therefore 
perturbative. However, this situation changes once |k| 
is large enough for the shift to reach the value u — Ij: 
above Kc the shift v splits into two imaginary solutions 
of the type 



v[k > Kc) = Ij ± «p(k) 



(86) 



Beyond this point, the renormalized angular momentum 
i'{k) cannot be expressed as a power series in terms of k, 
as a consequence of the non-analytical character of the 
split. Consequently, we expect the standard perturbative 
treatment to fail above this critical value Kc- The criti- 
cal reduced momenta for the uncoupled channels can be 
consulted in Table HI 



2. The Coupled Channel Case 

The (angular momentum) coupled channel case was 
worked out in detail by Birse [32|] . The extension is trivial 
and only requires the following definitions. In first place, 
the Schrodinger equation reads 



k.J 



T 2 



Ufe,,(r) = fc2^^fe,,(r), (87) 



where there are N angular momentum channels, Sj is the 
tensor matrix, and Lj is a diagonal matrix containing the 
value of the angular momenta 



\2 ^ diag(Zi(?i 1), . . . , In{In + 1)) • 



As in the previous case, the general solution of the 
Schrodinger equation is a linear combination of a the 
functions 

ufcj (^) = E ('^; fc) + Ph ^)] ' (89) 
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'Cj 




'Cj 




'-'Gj 










1.259 








1.328 


'Pi 


2.518 


'^Si-'Di 


1.367 




8.333 




4.314 




8.333 


'P2-'F2 


3.237 


''P2-''F2 


1.644 




1.089 




19.70 




7.913 




5.594 


^F:i-''Pi-''F3-''H-i 


2.615 



TABLE I. Reduced critical momenta Kc for tiie diflerent types of coupled channels appearing in the HH system {j < 3). 



where we sum over the possible values of the angular 
momenta. In this case ^ and r/ are iV-component vectors 
that can be expressed as sums of Bessel functions 



(90) 



m— — oo 

(91) 

where we have labelled the different solutions by the sub- 
script I ■ , which indicates that for k = the solutions will 
behave as free waves with angular momentum Ij. The 
bm{i^ij) coefficients are now vectors satisfying a certain 
recursive relation, see Ref. |32i for details. 

Of course, what is important are the analytical prop- 
erties of The functional dependence v = ^{k) is 
determined by finding the zeros of 



det iFj{iy, k)) =0, 



(92) 



where the characteristic function Fj{i/,K) is a N x N 
matrix. The equation above admits N solutions, one 
for each value of the angular momentum, which can be 
labelled by their value for k = 0, that is 



Vh = 0) = + - , 
with i ^ 1, . . . , N. We define Fj(i^, k) as follows 



(93) 



The analytical properties of the shift vi. (k) are anal- 
ogous to what happened in the uncoupled channel case: 
the shift begins its trajectory at i^j. = li + ^ for k = 
and moves slowly downwards. Once the shift reaches the 
value vi- — li at the critical value of the coupling k — k^, 
the shift splits into two complex conjugate solutions of 
the type I'liin) = U ±ipi-{K), signalling the breakdown of 
a long range perturbative expansion of the wave function. 

In general, the first shift to split into the complex plane 
is the one corresponding to the smallest angular momen- 
tum. For that reason we will simply study the critical 
value of K corresponding to the lowest angular momen- 
tum in the coupled channel. The results for the ^Cj, ^Cj 
and ^^^Cj families of coupled channels with j < 3 are 
listed in Table HI While the critical momenta for the ^Cj 
type of coupled channels was already determined by Birse 
in Ref. 32|, the results for the ^Cj and ^~^Cj coupled 



channels are new. 

In general, for the '^Cj channels we obtain values of Kc 
with are half the values obtained in Ref [s^] ■ The reason 
is simple: the tensor operator matrix elements in the ^Cj 
channels for heavy meson molecules is half of the value of 
the corresponding matrix in the two nucleon system. In 
general, the critical reduced momenta Kc grows rapidly 
with the angular momentum of the system, with the ex- 
ception of the "^Po partial wave for which a relatively low 
value of Kc is obtained. If we assume similar values of 
a3, the tensor length scale, the channel less likely to be 
perturbative is the ^ D2-^ 82-^ D2-^G2 partial wave, corre- 



sponding to a J 



PC _ 



P*P* meson molecule. 



¥j{y, k) = - — [Ro(i/) - Ro(-i^)] , (94) 
where fj {v) is a diagonal matrix defined by 

= Amg{v^ - {h + ^)^ ...,u^-[In + \f) ■ 

(95) 

In turn, Ro(i') can be expressed in terms of a recursive 
relation 



R„(j^) = [^j{v) - Sj R„+i(:^)] ^ 



(96) 



As in the uncoupled channel case, between 20 and 30 
iterations are more than enough to obtain a sufficiently 
accurate result. 



C. Finite Pion Mass Effects 

In the previous paragraphs we have calculated the crit- 
ical reduced momenta Kc for a pure \/r^ potential, cor- 
responding to the form of OPE for /i = 0. For /i > we 
expect the values of Kc to rise by a certain amount, as the 
finite pion mass e~^'^ will decrease the effective strength 
of the potential at large distances. Curiously, we will 
see that the dependence of Kc on /i is actually triggered 
by the far from perfect separation of scales in the chiral 
EFT, rather than from the finite pion mass alone. 

The argument is the following: we begin by solving 
the Schrodinger equation stepwise, that is, we divide the 
possible radii < r < 00 into several regions of a given 
size Ar. In particular, if we consider the region defined 
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by the condition 



VI. DISCUSSION AND CONCLUSIONS 



R 



Ar < r < R 



Ar 



(97) 



then, for fi Ar < 1 we can approximate the finite pion 
effects by the substitution 



0.3 



-fir 



0.3 



where we have ignored the 1 Ir and 1/r^ components of 
the tensor force for simphcitvrn. At this point it is worth 
noticing that the expUcit solutions of the Schrodinger 
equation we have written in the preceding subsections 
can also be applied in a specific subset of the real axis 
< r < oo. That is, if the 1/r^ potential is only valid in 
the region defined by Eq. ([97)) : 



Ar 

\R+—\-r) 



{r-\R-^\). (99) 



then, Eqs. dTS]), and ^ represent the full solution 
of the Schrodinger equation in this region. Moreover, the 
critical reduced momenta Hc calculated for < r < cx3 
are stiU vahd in i? - Ar/2 < r < i? + Ar/2. 

The consequence for a decaying e~^''/r^ potential is 
that, by making Ar sufficiently small, we can define a 
reduced critical momentum 



Kcifl.R) = KcC 



(100) 



which applies in the vicinity of R. Therefore, the real 
Kcip) setting the actual limits of perturbative tensor 
OPE is the minimum of KcifJ., R), which is in turn deter- 
mined by the smallest radius for which the OPE potential 
is valid. In chiral EFT we expect this radius to be propor- 
tional to the inverse of the breakdown scale, Rq oc I/Aq. 
Naively, we expect Rq to lie around 0.5 fm, but we cannot 
discard larger (or smaller) values. In particular, a recent 
analysis of the convergence of the two-nucleon potential 
in nuclear EFT [64*1 suggests values around 0.8 fm (and 
even higher). Of course, the extrapolation of this latter 
value to the heavy meson- antimeson potential is to be 
taken with a grain of salt, but we will nevertheless use 
the 0.8 fm figure only as an upper bound for Rq. For 
/i = m^r, the Rq = 0.5 — 0.8 fm window yields 



K(m7r) — Kci'm.Tr , Rq) — ( — a/S) Kc 



(101) 



where we have approximated e™^^° by \/2 (V3) for Rq = 
0.5 (0.8) fm. In the nucleon-nucleon case the a/2 factor 
approximately corresponds to the observed mismatch be- 
tween the critical momentum obtained by Birse for the 
^Si-^Di channel {pc = 66 MeV ) 13^ and the momentum 
at which the FMS calculations (Slf fails {pc ~ 100 MeV). 
However, a more recent formulation of nuclear EFT with 
perturbative OPE {6^ suggests a larger pc > 150 MeV, 
requiring Rq > l.Ofm at least. 



These components can be implicitly taken into account by a re- 
definition of 03. However, provided R is sufliciently small, they 
are suppressed by a factor of 11 R. 



The critical values Kc can be converted into critical 
momenta by dividing by the tensor length scale 03 



0.3 



(102) 



In the equation above, r — 21(1 + 1) — 3 is the eigen- 
value of the isospin operator ti ■ T2 and /ijjH the reduced 
mass of the heavy meson molecule under consideration. 
For the axial coupling constant g, we take g = 0.6 ± 0.1 
in the charm sector (D(*)D(*)) and g = 0.5 ± 0.1 in the 
bottom one (B(*)b(*)). For the charm mesons, the value 
of g is known from the D* — > Dtt and D* — > Dj de- 
cays [11,163. This yields g = 0.59 ± 0.01 ± 0.07, which 
we round to g = 0.6 ±0.1. For the bottom mesons, g is 
experimentally unknown, even though from heavy quark 
symmetry we should expect a value similar to the charm 
mesons one. There are however theoretical determina- 
tions, for example from lattice QCD [6^ - l7l] | or from the 
QCD Dyson- Schwinger equations [t^] to name just a few 
of them, usually ranging from 0.3 to 0.7 [zl]. The value 
we have chosen, g — 0.5±0.1, represent a compromise be- 
tween the previous determinations, but it should be kept 
in mind that the uncertainties for the bottom sector may 
be much larger. 

The critical momentum (including finite pion mass cor- 
rections) is defined as 



03 



(103) 



with Rq = 0.65 ± 0.15 fm (i.e. Rq = 0.5- 0.8 fm), and Kc 
to be taken from Table H] We can also translate Pcrit (m) 
into critical binding energies by 



-Scrit(A*) 



Pcrit(M) 

2^ 



(104) 



HH 



where we can make the approximation i?crit(™7r) — (2.5± 
0.5)i?crit(0). It is also interesting that the critical binding 
energy scales as l/g*. This dependence will generate 
remarkable uncertainties on the value of -Bcrit for bottom 
meson molecules. 

In the charm isovector sector, that is, the d(*^D'^*) 
molecules with I — 1 and t — I, the critical momenta 
are of the order of 0.8 — 1.1 GeV, similar to the chiral 
breakdown scale Aq ~ 0.5 — 1.0 GeV. In this case, cen- 
tral and tensor OPE are expected to be a perturbation in 
all the range of validity of the EFT. On the other extreme 
we have the bottom isoscalar sector, for which the critical 
momenta are of the order of nij^ (^ 200 MeV for p, = m-^). 
For these systems the onset of non-perturbative OPE 
start at around 4 MeV in the chiral limit, or around 
10 MeV if we take into account the finite pion mass ef- 
fects. It should be noted that the previous hmits could be 
more stringent if g is larger than the expected 0.5, as the 
critical binding depends on l/g*. However, the bottom- 
line is that in these heavy meson molecules OPE is ex- 
pected to be non-perturbative unless the bound states are 
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HH 


Pcrit(O) 


Bcrit(O) 


Pcrit(At) 




0±(l+±) 




DD* 


290t^g« 


42+46 


290t-» 


42lf, 


0+(0++) 

o-(i+-) 

0+(2++) 




D*D* 
D*D* 
D*D* 


270+]f 
220ir 


36+5? 
24l?« 


420tirc! 
44Q+200 

350ti«° 


qn+100 

f;Q+65 
OU_30 






BB* 






69ot«;; 


qo+130 

^"-50 


i-(o++) 

1+(1+-) 
l-(2++) 


^D2-''S2-''D2-''G2 


B*B* 
B*B* 
B*B* 


450t?i° 
360+?™ 


361?^ 
00+55 

•^8-20 
24+35 


cqo+390 
560^??" 


qn+130 



TABLE II. Critical value of the momenta for the tensor piece of OPE for the s-wave isoscalar DD* and D*D* and isovector 
BB* and B*B* systems. We compute the critical value both in the /i = limit and in the finite limit. In the second case, 
we usually have fi — m,r, with the exception of the DD* (BB*) case in which = 31 MeV (fi — 131 MeV), generating no 
significant deviation from /i = (/x = m-n)- The respective value of the critical momenta and maximum binding energies for 
the isovector DD* and D* D* molecules is three (nine) times larger than for their isoscalar counterparts. The contrary happens 
with the isoscalar BB* and B* B* molecules, for which the critical momenta are 1/3 of the isovector case and the maximum 
binding energies 1/9. 



genuinely shallow. The corresponding EFT treatment is 
the one employed in Ref. [28] for exploring prospective 
B*B states. 

The two interesting cases are the isovector charm and 
isoscalar bottom meson molecules: in these two systems 
the values of the critical momenta lie between m^r and Aq, 
defining a specific window in which OPE is perturbative. 
The critical momenta for the charm and bottom sectors 
can be found in Table HIl For the charm isoscalar the crit- 
ical binding energy lies around 25 — 40 MeV (60—90 MeV) 
ioi fj. = {fi = rriTr). In this regard, the D*D/DD* sys- 
tems are specially interesting in the sense that they natu- 
rally fulfill the condition ~ oE3- This means that the 
X(3872) is well within the energy range for which the LO 
EFT description consists only on contact range interac- 
tions. In particular, the critical momentum of 290 MeV 
implies that the charged D*"'"D^/D*~D+ component of 
the X(3872) can also be accounted for without the ex- 
plicit inclusion of pion exchanges at the lowest order. 
This represents a remarkable simplification that allows, 
for example, to reinterpret the explanation of Gamer- 
mann et al. [2^, [2^ for the branching ratio 

r [X(3872) ^ J/jjn+n-n^] 

r[X(3872)^ J/l^TT+TT-] ' ^ ' 

as a result within a specific EFT framework, or to extend 
X-EFT [lil to the charged channels of the AT (3872). 

In the isovector bottom sector, to which the recently 
discovered Zb(10610) and Zfc(10650) belong, the criti- 
cal energies range from 60 to 90 MeV, although they 
greatly depend on the axial coupling g, which value is 
not that well-known. This energy is in fact of the order 



Strictly speaking \fj,\ = 31 MeV m,r, which is small enough as 
to consider /i ~ as a pretty good approximation. 



of the splitting between the different channels (A^ = 
46MeV/92MeV), indicating that non-perturbative OPE 
is probably not needed until we reach binding energies 
such that coupled channel effects cannot be longer ig- 
nored. In this regard, the pionless EFT exploration of 
Mehen and Powell pLQ] probably provides a perfectly 
good LO description of the prospective HQSS partners 
of the Zfc's as far as their binding energies is less than 
40 — 50 MeV. Moreover, the inclusion of perturbative pi- 
ous in the subleading orders of the aforementioned EFT 
will be a direct extension of X-EFT [l^ to the bottom 
meson-antimeson system. We can call this theory Z-EFT. 

Of course, the previous results are to be held in proper 
context: they refer to the EFT framework rather than 
to a more fundamental or phenomenological level of de- 
scription of heavy meson molecules. In phenomenologi- 
cal approaches non-perturbative pion exchanges and cou- 
pled channel effects may be important in suggesting the 
possibility or explaining the location of heavy molecular 
states [6l|, [tI] , in particular if the states are shallow, in 
which case a fine tuning takes place between the different 
long and short range effects that are explicitly included. 
In contrast, in the EFT approach these fine tunings are 
absorbed in the counterterms, which are set to repro- 
duce the already known position of a certain state. This 
in turn stabilizes the calculations and reduces the impact 
of new physical effects we may subsequently incorporate 
in the EFT description as they are expected to be con- 
strained by power counting. The price to pay is the in- 
ability of the EFT framework to predict the positions of 
these states from first principles. 

To summarize, we have found that pion exchange ef- 
fects are perturbative in most types of heavy meson 
molecules over the expected range of applicability of 
heavy meson EFT. The exception is the isoscalar bottom 
meson molecules, for which OPE is particularly strong 
and presumably non-perturbative. To a lesser extent, in 
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the isoscalar charm sector non-perturbative OPE may be 
required if the molecular states are deep enough (around 
80MeV or more). At lowest order, below the critical 
binding energies we have computed, the EFT description 
is a contact-range theory. Surprisingly, this description 
holds for momenta much larger than naively expected. 
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Appendix A: Derivation of The Leading Order 
Potential in Heavy Hadron Chiral Perturbation 
Theory 

In this appendix, we derive the LO potential between a 
heavy mesons and anti-meson within heavy hadron chi- 
ral perturbation theory (HHChPT) 01 . The LO po- 
tential can be decomposed into a contact-range and a 
finite-range piece. The later can be identified with the 
well-known one pion exchange (OPE) potential. We also 
compute the complete partial wave projection of the LO 
potential. 



depending on whether we are considering the pseu- 
doscalar or vector meson case. 

The HHChPT Lagrangian at lowest order {Q°) can be 
decomposed as the sum of two different contributions: 



(A7) 



where the first term represents a 4-meson interaction ver- 
tex and the second the meson-pion vertex. The contact 
interaction term reads nm 



C 



12 = i^OaTr [hQrQ^^] Tr [H'^H'ir 
+ Dob Tr [H«H«7^75] Tr [h^ 11^^15 
+ Eoa Tr [i7«T,H«7^] Tr [i/'3T,5'37^ 
-f £^06 Tr [i?«r,77«7^75] Tr [i/'35'3r,7^75 



(AS) 



where are the isospin matrices (i.e. the Pauli matri- 
ces). That is, in the previous Lagrangian D^a and Doh 
are isospin independent, while E^a and Eqi, are isospin 
dependent. In general we will be only considering a spe- 
cific isospin channel, in which case we will generically 
define 



Coa — +Ti ■ T2 Eoa , 

Cob = Dob + n • T2 Eob ■ 



(A9) 
(AlO) 



For the Lagrangian describing the pion-meson vertex we 
have El 



1. The Effective Lagrangian at Leading Order 

We define the heavy meson (antimeson) fields in terms 
of the pseudoscalar and vector meson (antimeson) fields 
in the following way. 

1 + f 



HQ = joH^ \o ■ 



(Al) 
(A2) 

(A3) 
(A4) 



where 7^ and 75 are the Dirac gamma matrices, v is 

the velocity parameter and f = w^7^. The field P^*^-* 

(Pi*'^^) destroys a heavy meson (antimeson) of velocity v 
and polarization /j, for the vector meson. The subscript „ 
also indicates that the normalization of the one-particle 
states is given by 

(P„(fc)|P„.(fc')> = ^v"5,,, (2^)3<5(3)(fc- fc') , (A5) 
(P7(^)|P,r (fc')) = 2i;° 5,,, (2^)3 <5(3)(^ _ k') s^, , 

(A6) 



C 



(0) 



I — Tr 
2 



where ^ is defined as 









M = 




V5 



V2 



(All) 

(A12) 
(A13) 



with vr the pion fields and /^r = 132 McV. The isospin 
indices convention for the heavy meson and antimeson 
fields is 



P 



,P 



1",+ 
(*m) 



and Pi*'^) 



-P, 



u,0 



(A14) 



This convention explains the minus sign in the second 
line of Eq. (jAlip . That is, we have performed a G-parity 
transformation to obtain the chiral Lagrangian for the 
heavy anti-meson fields from the heavy meson one. 
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The Non-Relativistic Normalization 



3. The Contact Range Potential 



In the heavy quark limit, toq oo, the static poten- 
tial between two heavy mesons is a well-defined object. 
Relativistic effects are suppressed, and the two heavy me- 
son system can be effectively described in terms of non- 
relativistic quantum mechanics. For making the non- 
relativistic transition we specify the velocity parameter 
to be I) = (1,0) and employ a new normalization of the 
heavy meson fields defined by 



(A15) 



where we have substituted the Greek index ji by the latin 
index j, as the /x = polarization component is com- 
pletely irrelevant for t; = (1, 0). In this heavy meson field 
normalization, the one-particle states have the usual non- 
relativistic normalization 

{P{k)\P{k')) = (27r)3,5(3)(fc-fc'), (A16) 
{P*\k)\P*^{k')) - {2nf5'^''\k-k')5,, , (A17) 

which is more convenient for the definition of the quan- 
tum mechanical potential. 



In the non-relativistic normalization, after expanding 
the H*^*^^ and H'^'^) fields, the contact range Lagrangian 
can be rewritten as 

Ah (a) = -Coa {P^P + P* ^P*) {PP^ + P*P* ^) 

+Cah (P* P^ +PP*^)iP* P^ +PP*^) 
-iCob\{P* pt +pp*t-)(p*t X p*) 

- (P*t X P*){P* Pt +p_p*t 

+Cob {P* ^ X P*) (P* t X P*) . (A18) 

In the expressions above, the polarization of the vector 
mesons has been made implicit to simplify the notation. 
The potential corresponding to the Lagrangian above can 
mix different particle channels. Thus, it is convenient to 
write the potential in the coupled channel basis 

Shh = {l^'^), 1^*^), l^'^*), 1^*^*)} , (A19) 
for which we have (in momentum space) 



Vc{q) = C, 



Oa 



/l 

ei* -el 

€2* ■ el 

\0 ei*-e{e2' 



06 



-er*-eJ/*\ 

-el • €2* +Si ■ €2'* 

— ei'* • €2 — eT'* • S2 

-el • €2 +Si ■ €2 -el • ^2 • S2 ) 



(A20) 



In the equation above, t\{2)l^v(2') represents the polar- 
ization of the incoming particles 1(2) or the outgoing par- 
ticles l'(2'). In turn, the matrix elements of the spin-1 
operator S\{2) are equivalent to the vector product of the 
polarization wave functions, that is 



i(lA'|5|lA) = (ev*xeA), 



(A21) 



where |1A^'^) is the vector corresponding to the polariza- 



tion wave function e_j^(/) . 



4. The Finite Range OPE Potential 

For obtaining the finite range piece of the potential, 
we start by considering the HHtt vertex in the non- 
relativistic normalization, for which the Lagrangian reads 



■'■kHH 



[{P*P^ +P*^P) + i (P*t X P*)] f • a? 



+ ^i— [(P*Pt + P* tp) - i (P*t xP*)]T-dTT, 

(A22) 



where the polarization indices are implicit. If we define 
the non-relativistic amplitude as follows 

A{H ^ Htt'') = i {H\£„HH\H7r'') , (A23) 



we find 



A{P Ptt") = , (A24) 
A{P*' ^ Ptt") = ^^ei-q, (A25) 

/tt v2 

AiP*' P*%") = -ij-^S-q, (A26) 

where q is the momentum of the outgoing pion, a the 
isospin index of the pion and e^ the polarization wave 
function of the heavy vector meson P**. The amplitudes 
with heavy anti-mesons can be obtained either from the 
Lagrangian or from a G-parity transformation of the ex- 
pressions above, yielding 

A{P Ptt") = , (A27) 
A{P» ^ Ptt'^) = f ^ el • q, (A28) 

A{P*' P*%") =ij-^S-q, (A29) 
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where only the P*P*7r vertex changes sign. 



From the non-relativistic amphtudes above, the OPE 
potential in the heavy quark limit is given by 



Tl • T2 







'Ey* ■ qt2- q 
f 5i ■ qe^2- q 





-el • qei* ■ q 


-el • (f 5*2 • q 



-ev ■qe2' ■ q\ 
+Si ■ qei* ■ q 
-ey* ■ qS2- q 
+Si ■ qS2- q 



(A30) 



where we have already particularized for the heavy 
meson-antimeson case. If we consider the effect of the 
splitting between the masses of the PP, PP*/P*P, and 
P*P* systems, then we should substitute the pion prop- 
agator by 



(A31) 



in the non-diagonal terms of the potential matrix, where 
the effective pion mass is given by /^^ = — Ag, with 
Aq the energy splitting. This corrections are fundamen- 
tal in the charm sector, where Aq is of the order of 
m^. A further problem arises in the PP — >■ P*P* and 
PP*/P*P — > p*p* transitions, where either the initial or 
final heavy meson states will be off-shell. In this case it is 
not entirely clear whether the static approximation can 
be employed to define a transition potential between the 
different heavy meson states in the charm sector. Luck- 
ily this problem is not present in the heavy quark limit, 

I 



I 

in which the possible off-shell effects are strongly sup- 
pressed. 

We can calculate the coordinate space representation 
of the OPE potential by Fourier transforming the mo- 
mentum space potential 



from which we obtain 



^ V^°^ - in) e''^-^ 
(27r)3 '^^'.hhW'' 



(A32) 



V, 



(0) 



(f)=f2.fiCi2(f)-^5(f) 



F.HH^' ' —i-^v I Qj2 



T2 ■ Tl 



Ci2{f)Wc{r) 



o{—), 



Si2{f)WT{r) 

(A33) 



with Wc{r) and Wrir) defined in Eqs. ^ and ((571) . 
The central and tensor operators Ci2(r) and Si2(r) are 
4x4 matrices in the particle channel basis yBjjg. Their 
explicit representation is 



Ci2(r) 



Sl2(f) 



f 

-ei-e2'* 

-eT'*-e2 

\-el • 62 +5*1 • e2 -el • 5'2 

/ 





-ei' • £2' \ 
-ei'* ■ S2 

+ Si-S2 J 

-5'i2(er'*,el'*,f)\ 
-S'i2(el,e2'*,f) +Si2{Si,e2'* ,f) 



~Si2{ey* ,£"2,f) ~Si2{<^i'*,S2,f] 

\-S'i2(el,e^,f) -fS'i2('5'i,e2,f) - S 12(^1, S2,r) +Si2{Si, S2,r) J 



(A34) 



(A35) 



where S12 is defined as 

<S'i2(ai, ^2, ^) = 3di ■ fb2 ■ r - Si ■ b2 

5. Partial Wave Projection 



(A36) 



We can simplify the representation of the OPE poten- 
tial by projecting into Hfl states with well-defined J^*^ 



r 



quantum numbers. This tasks is greatly simplified in co- 
ordinate space, for which only the operators Ci2{f) and 
Si2(r) are to be projected. In the PP case, we consider 
the partial waves 



\PP{ljm)) = SjlYlrr^f) 



(A37) 



where Yi„i{f) is a spherical harmonic. The parity and 
C-parity of this partial wave are P = C = (—1)'. The 



19 



PP* / P*P cases are a bit more involved owing to the po- 
larization of the P*/P* vector meson. The partial waves 
with well-defined orbital and total angular momentum I 
and j are defined as 



\P*P{ljm)) 



E 

mi,n 



Yirr,,{f)\lu){lmilu\jm), (A38) 
Ylrr,,{r)\ln){lmll^i\Jm) , (A39) 



where |l/u), are the polarization vectors of particle 1 
and 2 and {Ii'mil2m2\jm) a Clcbsch-Gordan coefficient. 
The polarization vectors are to be understood as giving 
matrix elements of the type 



(V|ep*-ei|l/x) 
(lzv'|e2'*-ei|l/i) 



(A40) 
(A41) 



and so on, where we have added primas to denote final 
states and with the unit vector in the spherical basis. 
A problem with the two partial waves we have written 
above is that they do not have a well-defined C-parity. 
Thus, we define the linear combination 

\PP*{rj){ljm)) = i= [\PP*{ljm))-ri\P*P{ljm))] , 

(A42) 

with intrinsic C-parity 77 and total C-parity C = (—1)' tj. 
Finally, the P*P* partial wave with intrinsic, orbital and 
total angular momentum s, I and j is given by 

\P*P*{sljm))= ^ Yirni{f)\s'ms){l'mis'ms\3m) , 

mi,ms 

(A43) 

where the spin state \sms) is 

\sms) = ^ |l/i)|l!^)(l/ul!^|sms) . (A44) 



The C-parity of these states is C = (—1)'+*. Alterna- 
tively we can also use the spectroscopic notation "^^^^Ij 
instead of {s)lj for denoting the partial waves. 

From the previous definitions we can compute the ma- 
trix elements of the central and tensor operators as 

((s')?'/m'|6i2|(s)/jm) = 

d^f {{s')l'fm'\6i2{f)\{s)ljm) = 



(A45) 



where Oi2('r) either represents Ci2(r) or Si2(r). The 
particle channels have been made implicit owing to the 
matrix notation. The central and tensor operators pre- 
serve the J^'^ quantum numbers, which for the particular 
case of the tensor operator means that the orbital angular 
momentum can only change by an even number of units. 
In addition, if we only consider the P*P* particle chan- 
nel, the intrinsic angular momentum is also subjected to 



the restriction of even |s — s'|. However, terms mixing the 
P*P* and pp*/p*p channels can change s by one unit. 

Thus the most compa,ct way to write the matrix ele- 
ments of the Ci2(r) or Si2(r) operators is to consider a 
particle/partial wave basis with well-defined J^*-^. Ex- 
amples are the following 

B^n{0++) = {\'SoiPP)),\'S„iP*P*)), 

~MP*P*))}, (A46) 



B^ji{l+-) = [fSi{PP* + P^P)) 



\^Di{PP* +P*P)), 

\'SiiP*P*)),\'D,iP*P*))},iAA7) 

Bhu{2++) = {\'D2{PP)), \'D2{PP* - P*P)), 

|l£)2(P*P*)),|^52(P*P*)), 

fD2{P*P*)), fF2{P*P*))} , (A48) 
for which we obtain the central matrices 



Ci2(l + -) 



Ci2(2++) 



and the tensor matrices 





73 0\ 




V3 -2 , 




1/ 


/-I 


2 \ 





-10 2 


2 


0-10 


Vo 


2 -1/ 


/o 


V3 o\ 





1 


V3 


-2 





10 





1 


\o 


1/ 



(A49) 



(A50) 



(A51) 



&12'" 



(0++) 



Sl2(l+-) 



^ -Vg 

-V2 
y-V6 -V2 -2 

/ V2 V2\ 

\^ -1 \/2 -1 

\/2 \/2 

\V2 -1 V2 -ij 



(A52) 



(A53) 
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Other J^'-^ possibilities are to be computed as these three 
examples. As can be seen, the resulting coupled channel 
structure can be quite complex: in certain cases we can 
have up to six coupled channels. The simplifying feature 
is that particle coupled channel effects are suppressed by 
two orders in the chiral expansion of the potential, which 



means that in low order calculations we can concentrate 
in a specific particle channel and ignore the others. In 
this case, we only need to consider the diagonal terms of 
the central operator and, in what regard the tensor op- 
erator, we end up with the matrices shown from Eq. (|5ip 
to (|61l) with at most four coupled channels. 
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